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Abstract 



The aim of this paper is to derive explícit formulae for the Riemannian Bures 
mètric g on the manifold T> of (finite dimensional) nondegenerate density ma- 
trices g. The computation of the Bures mètric using these equations does not 
require any diagonalization procedure. The first equations we give are, essen- 
tially, of the form g = Y ajjTr dg g l ~ l dg g 3 , where üíj is a matrix of invari- 
ants of g. A further formula, g = ^ Cij dpi (g> dpj + Yl bijTr dg g 1-1 dg g- 1 " 1 , is 
adapted to the local orthogonal decomposition T> ~ W 1 x u ^/t" at genèric 
points. 
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I. INTRODUCTION 

Let D be the manifold of all nondegenerate, positive, Hermitian nxn- matrices. The 
tangent space T e T> consists of all Hermitian nxn-matrices and the Riemannian Bures mètric 
on V is defined by, [p]], 

g e (Y',Y) = ^TrY'X, Y, Y' E T g V , (1) 

where X is the (unique) solution of qX + Xg = Y. The submanifold of trace one matrices 
is the space of so called completely entangled mixed states of a finite dimensional quantum 
system. This mètric appears quite naturally on the background of purifications of mixed 
states and is used in quantum information theory to describe the statistical distance of mixed 
states. It is an extremal monotone mètric and seems to be quite distinguished for several 
mathematical and physical reasons, see e. g. @. It should be mentioned, that ([[]) defines 
also a mètric on manifolds of stable range densities, but we will deal with the maximal range, 
only. 

Several formulae and approaches for computing the Bures mètric were given for low 
dimensions, e. g. 0,^-0] • li \a) , a = 1, . . . ,n, are eigenvectors of g with eigenvalues X a , then 
a simple calculation shows that ([[]) yields, |§ , 

a,/3 

but such a formula we will not accept as explicit, because it needs the knowledge of eigen- 
values. 

The aim of this note is to provide several equations for computing the Bures mètric 
in any finite dimension using matrix products, determinants and traces, only. Section II 
provides expressions for the Bures mètric based, essentially, on results of the theory of matrix 
equations. In Section III we give an equation adapted to the local isometric decomposition 
of the manifold T>. 

Notations: The following quantities depend on a positive, Hermitian n x n - matrix g. In 
order to simplify the notation the dependence on g will be suppressed. By Ai < . . . < À n we 
denote the eigenvalues of g and by A the diagonal matrix of the Àj-s. Moreover, V will be 
the Vandermonde matrix (A^ 1 ) . Operators acting on matrices are denoted by bold letters, 
especially, if not indicated otherwise, L and R denote the left and right multiplication by g. 
The Bures mètric now takes the form 

We set 

X{t) := det (íl - q) = t n + M n_1 + • • • + K, 

ko := 1 and ki := for k > n or k < 0. Hence, := (— l) 1 ki is the elementary invariant 
of degree i and (— l) n x(t) is the characteristic polynomial of g. Moreover, let Pi := Tr g % so 
that the differential dpi applied to a tangent vector X yields dpi(X) = iTiXg l ~ l . Finally, 
we will several times make use of the matrix 



2 



P :-- 



Pl P2 
P2 P3 



Pn 

Pn+1 



p n p n+1 . . . p 2n -l 

For instance we have the 

Criterion: A Hermitian nxn-matrix g > is a genèric point of T> iffdetP ^ 0. 
Indeed, P = V T AV. Therefore, det P = det(A) det(V) 2 = det(g) U t<j ( X i ~ X j) 2 - 



(3) 



□ 



II. GENERAL FORMULAE 



In order to calculate g e (Y', Y) one has to solve the matrix equation 

gX + Xg = Y . (4) 

Then 2g g (Y',Y) = TrY'X. The matrix (resp. operator) equation EX — XF = Y was 
intensively studied (for a review see 0). A bàsic item is that it has a unique solution X if E 
and F have disjoint spectra (Sylvester, Rosenblum). In our matrix case, E = g = —F, this 
is fulfilled because g is positive. The uniqueness is also clear, not appealing to this theorem, 
since we may suppose w.l.o.g. g to be diagonal with eigenvalues A« and equation (g) becomes 
(Aj + Xj)Xij = Yij] i,j = l,...n. A further, simple but nice, tool in this theory is the use of 
similarities of block matrices, in our case e. g. 



" 1 


-x~ 




-Q 


Y ' 




1 


X' 




-Q 


o" 





1 







Q _ 







1 







Q _ 



(5) 



Applying the polynomial x as an operator function to both sides leads in the upper right 
box to the identity 

X (-g)X + M = 0, (6) 
where M is the upper right box of x applied to the above middle matrix containing Y; 

n i—í 

But x( — o) is invertible. This can be seen as follows. The characteristic polynomial of — g 
equals (— l) n x(— t) and the positivity of g implies that x(í) and x(— t) have no common 
divisors. Hence, by the Euclidean algorithm, there exists two polynomials p, q such that 
p(t)x(t) + q(t)x(—t) = 1, and inserting g gives q(g)x(~ o) — 1- Therefore, the solution X is 
given by 

n i — í 

x = -x(-q)- 1 E k ^ E(-^ iy ^' _1 » ( 8a ) 

i=í 3=0 

or, in a more compact form, 
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1 X 





x(-qY 





A 



■q y 

Q 



(8b) 



The first explícit formula we get for the Bures mètric is: 
Proposition 1: 



g(Y',Y) :=-^Tr 




Y' 



x(-q)- 1 o 





X 



-Q y 

Q 



(9) 



□ 



The inverse of x(~ o) is again a polynomial expression in g. Therefore, rewriting (8a) using 
the Cayley-Hamilton theorem, X can be transformed to the form Yli<ij< n a ijQ t ~ 1 y Q^ 1 with 
exponents less than n and coeficients being invariants of g. That means, 



1 



L + R 



E 

l<i,j'<n 

(Id, L, . . . , L"- 1 ) A (Id, R, . . . , R"" 1 ) 1, , A := (a 



h3, 



(10a) 
(10b) 



(Of course, this can be seen also directly. Since L + R acts on a finite dimensional space, 
its inverse is a polynomial in L + R, which again can be reduced to the above form using 
x(L) = L x (e) = anc i similar for R.) The solution of (f|) given by Smith in |§ is, essentially, 
of the above form YIkí j<n o·íjQ t ~ 1 YQ^~ 1 · The formulae for the coefficients in terms of 
invariants of g one reads off from || will be given at the end of this section. 

The representation ( |1~ÜD is unique provided g is genèric. Indeed, if there were coefficients 
such that X^kí j<n a ú'I jl_1 R- ?_1 = then applying this operator to all vectors of a common 
eigenbasis of L and R would result in V T AV = 0. But the Vandermonde matrix V is 
nondegenerate for a genèric g and we would conclude = 0. Moreover, in the genèric case 
the matrix A is necessarily symmetric. 

To get a compact expression for the coefficients in ([H]) we define a nxn-matrix K by 



K :-- 










-kr. 



1 







1 












-hi 



(11) 



K carries out the reduction of powers of g by x(o) — 0- Indeed we have 







' 1 1 


g 




g 




= K 




. g n . 







and similar for the reduction of powers of L and R. Thus, multiplying ( |10b|) by L + R leads 
to 



K L A + AK = C. 



(12) 
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where 



C 



1 ... 
... 



... 



(13) 



represents the identity operator Id = L°R°. Note that K has the same characteristic 
polynomial as g. Now we may proceed as above to find A. We replace in (|5[) —g, g, X and 
Y by — K T , K, A and C and apply Instead of (^) and (0) we obtain 



A = - X {-K 



T\-l 



N. 



(14) 



where 



i-1 



N = X 





c' 





K 


ün— 1 




k n -2 




IY' 1 






12 



3=0 



i=l 



fel 1 

-1 







-l) m fe 



n+Y—i—j 



J i j'=l 



(15) 



To see the last equation note that K Tk CK^ has only a 1 in the (fe + 1, j + 1) - position; K 
moves the 1 coming from C to the right and K T moves it down. Hence we get 

Proposition 2: The Bures mètric equals 



1 n 

g = - ^ a ij Tr d @ É? ? 1 d Q 



(16) 



where (a^) = — x(— ií T ) _1 A^ , fí and AT gróen 6y fcLdJ and (fUL). 
The probably "most explicit" form of the coefficients is given by 
Proposition 3 (R. A. Smith): 



□ 



n—i n—j 



2 detH 



^^(-1)^$ 



■ i + j + r + s , 



(17) 



r=0 s=0 



where 





"fei 


fes • 


• fe2n-l 


H = 


feo 


fe 2 • 


• fe2n-2 








. 


. 


fen 



"2,3 -i Jí,7 = l 



(18) 



and $(m) = if m is not an integer, and otherwise $(m) is the cofactor in det H of the 
element in the first row and m-th column of H . □ 
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Remark: The determinant of H is not equal zero, raore precisely: 

det H = (-I)t fi Xi l[(X t + Xj) + 0. 

i i<j 

n(n + l) 

Indeed, changing the order of rows and columns we get det H = (— 1) 2 det [e n+ i_2i+j] . 
But the last determinant is just the (symmetric) Schur function of the eigenvalues of g 
(comp. 0, 1.3) related to the partition (n, n — 1, . . . , 1) = (1, . . . , 1) + (n — 1, . . . , 0) leading 
to the above product. □ 



III. A FORMULA ADAPTED TO V w R n x u (") / T « 

Every g e T> can be diagonalized with a suitable unitary u; 

g = ufi 2 u* , /i = diag(/ii, ...,//„), G R+ 
and we have = 2u fidfiu* + u [u* du , /i 2 ] w*. The inverse problem now splits into 

- — 1 —— (dg) = u ( //~ x d// + ^— — ([w*dw,/i 2 ]) J w* 

and we get for the mètric 

g e = Tr d/x d/i + ^Tr [u*du,fi 2 ] 1 ([u*dw,/i 2 ]) . (19) 

Therefore, in a neighbourhood of a genèric point the Riemannian manifold T> locally looks 
like W 1 x u ^"-'/t") where W 1 is equipped with the Standard mètric and the mètric on the 
homogeneous space u (™)/t™ depends on the first parameter. The tangent space at g splits 
into 

T e V = T , e + r^ L· , (20) 

where T^ is the subspace of Hermitian matrices commuting with g. Its orthogonal comple- 
ment (w. r. to the Bures mètric) is the space of all [a, g], a - antihermitian. If g is diagonal 
then (|20|) is the decomposition into diagonal and off-diagonal Hermitian matrices. 

From now on let g be a genèric point. By P and = Id — P we denote the (orthogonal) 
projectors onto the subspaces in (|20|) . 

Lemma: 

n 

P(Y)= (P^.Triy- 1 (21a) 

n 

= E ( p %^ '• ( 21b ) 
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Proof: To show ([21 a| ) we use that for a genèric g the powers g, g 2 , . . . , g n form a basis of the 
vector space of all Hermitian matrices commuting with g. Moreover, 4P i .,· = g e {g 1 ., g') and 
4 Tr Y gi~ l = g s (Y, g>). Having this in mind ( |21a| ) is just the usual formula for the orthogonal 
projection onto a subspace with a given basis. To see ( [21b|) let Y a p, a,/3 = 1, . . . , n be a 
common eigenbasis of L and R, gY a p = X a Y a p, Y a p g = \pY a p (Y a p - may not be Hermitian). 
Then the complex span of T^ resp. T^ is generated by all Y a p with a = (3 resp. a ^ (3. For 
y — Yap the right hand si de of ( |21b| ) yields r\ a (iY a p, where 



V - the Vandermonde matrix of eigenvalues of g. But P = V T A V implies r] a p = 5 a p. □ 
From (|2l|) we now get 



J— p(y) = i y t l (p- 1 ) ..^(Y), 

^P^) = ^{a (i -i(p-') i .} í ,^-' 

i,.7=l 



where we used Trj'F^ 1 = dpj(Y) and l^r;(F) = |í? l Y for F G T". The matrix (a^) is 
given by Proposition 2 or 3. Inserting these equations into 

g = ^Tr (p(dg) J -^P(dg)+ dg^-^P^dg) 

yields 

Proposition 4: T/ie decomposition g = + g\ j_ o/ í/ie Bures mètric is given by 

•=ltfn^ + J £ (2«« Tïde^ 1 ^. (22) 

□ 



IV. EXAMPLES 



Propositions 1-4 involve elementary and power invariants, which can be expressed by 
each other (comp. |§). Rewriting the identity 



mk m + Y^rkm-r = , 771=1,2,... 



r=l 

as a system of linear equations for the e^-s resp. the pj-s one gets the relations 

i [e! 10. 

2e 2 ei 1 . 



det 



Pi 1 

P2 Pi 2 



Pi-1 Pi-2 Pi-3 • • • «"I 
Pi Pi_! Pi_2 ... Pl 



Pi = det 



(i-l)ej_i e^ 2 ei- 3 
ící ej_i ei_ 2 








1 

ei 
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This identity also allows for expressing the power invariants p m , m > n, by invariants of 
degree less or equal n. Especially the i + 1 - row of our matrix P equals (pj+i, . . . ,p n +i) = 
(p 1 ,...,p n )K i . 

The number of terms in Propositions 1-4 rapidly increases with the dimension n. Thus 
we give only certain expressions for n = 2, 3 in terms of power resp. elementary invariants. 
Concerning Proposition 1 we get 

1 ( TrF V + e lQ + e 2 l)- 1 (gY -Yg + e x Y) for n = 2 

g(Y',Y) := - \ TiY'(g 3 + e 1 g 2 + e 2 g + e 3 l)- 1 x 

[ (Q 2 Y ~ QYQ + YQ 2 + ei{gY - Y g) + e 2 Y) for n = 3 

The following terms appear in Propositions 2-4: 
n = 2: 



1 



9U 



1 

4e 2 (ef -4e 2 ) 



/ dpi dp_2 \ 
V 1 ' 2 / 



( dei , de 2 ) 



-P2 

-P2 Pl 

eie 2 -2e 2 
-2e 2 ei 



dpi 

d£2 
2 



dei 
de 2 



n = 3: 



2A - P" 1 



ei(eí -4e 2 ) 



2eie 2 

-2e 2 e 1 
ei -2 



ef + e 2 -ei 
-ei 1 



Pi(2p 2 -pi) 



Vi - P\ Pi 
Pi -2 



det P = -pjj + 2p 2 p 3 p4 - Pip\ - pIpò + PiP3Pò 



detP 



- P| - P2P5 P2P4 - P3 



* 
* 



P1P5 - P3 P2P3 - pm 
* pm - p\ 



A = 



1 



2e 3 (eie 2 - e 3 ) 



eie 2 + eie 3 - e 2 e 3 -ef e 2 eie 2 - e 3 



* 



ef + e 3 
* 



ei 
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